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UNCLASSIEIED 


Intrabeam  Scattering  and  Halo  Formation  in 
Intense  Ion  Beams 


N.  Uhlmann*  and  G.  Zwicknagel* 


* Institut  fiir  Theoretische  Physik,  Universilat  Erlangen-Nurnberg, 
Staudtstr.  7,  D - 91058  Erlangen,  Germany 


Abstract.  An  essential  problem  for  the  operation  of  high  current  linacs  is  the  control  of  beam 
losses  due  to  halo  formation.  Here  we  concentrate  on  an  interplay  between  intrabeam  scattering 
(IBS)  and  the  incidence  of  particles  which  are  driven  to  high  amplitudes  by  resonances  with  the 
nonlinear  space  charge  fields  of  a mismatched  beam.  For  these  studies,  we  developed  a particle- 
core-molecular-dynamics  (PCMD)  method  which  suitably  joins  the  time  evolution  in  framework 
of  the  envelope  equations  and  a calculation  of  IBS  between  pseudo-particles  with  a renormalized 
charge.  With  this  technique  we  investigated  continuous  KV-beams  in  a FODO  channel.  Our  first 
results  strongly  suggest,  that  IBS  in  fact  supports  the  halo  formation  in  a mismatched  beam. 


INTRODUCTION 

One  of  the  most  critical  aspects  for  a successful  operation  of  a high  current  linac  is  the 
control  of  beam  losses,  since  only  a very  small  fraction  of  particles  can  be  tolerated  to 
collide  with  the  accelerator  structures.  Such  beam  losses  are  associated  with  the  forma- 
tion of  a halo  of  particles  significantly  outside  the  beam  core.  Substantial  theoretical 
effort  is  presently  addressed  to  the  understanding  of  halo  formation,  e.g.  to  the  stability 
of  collective  excitation  modes  which  has  been  studied  in  the  framework  of  the  nonlin- 
ear Vlasov-Maxwell  equations  by  both  analytical  and  numerical  approaches  [1,  2,  3,  4]. 
Also  the  resonances  between  nonlinear  space  charge  forces  outside  the  core  of  a mis- 
matched beam  and  the  betatron  oscillations  have  been  identified  as  a source  for  a halo. 
This  has  been  studied  mainly  by  a frequency  map  analysis  of  the  single  particle  motion 
in  framework  of  a particle  core  (PC)  model  [3],  Hard  collisions  between  the  beam  par- 
ticles, which  are  neither  included  in  the  Vlasov-Maxwell  equation  nor  the  particle  core 
approach,  could  also  scatter  particles  into  a halo.  But  for  the  typical  beam  parameters  of 
high  power  linacs,  it  is  very  unlikely  that  this  intrabeam  scattering  (IBS)  is  sufficiently 
effective  to  directly  drive  a beam  halo  [5].  IBS  may,  however,  play  a non-negligible  role 
in  amplifying  the  other  mechanisms.  We  here  focus  on  an  interplay  between  intrabeam 
scattering  (IBS)  and  the  incidence  of  particles  which  are  pushed  to  high  amplitudes  by 
nonlinear  resonances  driven  by  the  space  charge  fields  of  a mismatched  beam.  As  a fully 
microscopic  numerical  treatment  including  all  mutual  Coulomb  interactions  of  the  beam 
particles,  and  thus  all  space  charge  phenomena  and  IBS,  considerably  exceeds  the  ca- 
pacity even  of  modern  computers,  we  developed  a alternative  method.  These  so-called 
particle-core-molecular-dynamics  (PCMD)  simulations  suitably  join  the  time  evolution 
of  the  beam  core  in  framework  of  the  envelope  equations  with  a microscopic  calculation 
of  the  Coulomb  interactions  between  pseudo-particles  with  a renormalized  charge. 
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THE  PARTICLE-CORE-MOLECULAR-DYNAMICS  METHOD 


The  equations  of  motion  for  the  position  ra  and  the  velocity  Va  of  a particle  with  mass  m 
and  charge  q in  an  accelerator  structure  or  a storage  ring  can  be  written  in  a local  center 
of  mass  frame  of  the  beam  as 

% = Va,  ^ = K(r«,/)  - -Vtl)(rc(,/).  (1) 

dr  dr  m 

Here  K denotes  all  the  forces  on  the  particle  due  to  the  beam  manipulating  elements  of 
the  structure  and  ^ represents  the  potential  generated  by  all  the  other  charged  particles 
in  the  beam.  In  case  of  an  ideal  FODO-channel  extending  along  the  z-dircction,  the 
external  forces  are  linear  and  K(r.r)  = {Ki{t)x,K2{t)y,0),  where  the  A',(r)  only  depend 
on  timer  or  likewise  on  the  position  .v  along  the  channel.  In  general  a numerical  treatment 
of  Eq.  (1)  is  unavoidable.  The  most  complete  approach  are  molecular  dynamics  (MD) 
simulations  which  take  completely  into  account  all  the  mutual  Coulomb  interactions 
between  the  beam  particles  and  thus  include  IBS  without  any  approximation.  They  have 
been  applied  successfully  for  low  current  beams  and  studies  of  e.g.  beam  crystallization 
and  related  issues  [6,  7,  8,  9,  10).  Because  the  numerical  effort  scales  like  N~  the 
computational  effort  is,  however,  prohibitive  for  high  current  beams  which  require  at 
least  ~ 10^. 

If  one  is,  on  the  other  hand,  interested  e.g.  in  the  motion  of  a few  particles  in  tlie 
nonlinear  space  charge  fields  outside  the  beam  core,  one  can  assume  a given  space 
charge  field  due  to  the  large  number  of  core  particles  with  regular  motion  and  ne- 
glect the  charges  of  the  actual  test-particles.  This  approach  becomes  easily  applica- 
ble for  the  Kapchinsky-Vladimirsky  (KV)  distribution,  where  a fairly  simple  .self- 
consistent  solution  exists.  The  continuous  KV-beam  has  the  uniform  charge  density 
p(s)  = 0(  1 -x^/a-{s)-y-/b-{s))/nci{s)h{x)  with  envelopes  «(.v)  and  b{x)  which  vary 

with  the  position  in  the  FODO  channel  .v  ■ypef  according  to 

d^ci  E 

— -y  -I-  /f,(.v)a  - — = 0 and  x y,a  b . (2) 

ds^  ^ ' a + h cd  ^ ' 

In  Eq.  (2)  the  constants  e,  arc  the  emittances,  E = ((//)/ (Tteowte^Y^P^)  is  the  space 
charge  parameter  for  a beam  current  /,  X-  = //ctP  is  the  line  density  and  P,y  are  the 
relativistic  factors.  Here  the  potential  generated  by  the  space  charge  of  the  beam 

can  be  obtained  analytically  from  the  known  p(/)  via  the  Poisson  equation.  Combining 
now  the  time  evolution  of  the  core,  given  by  a numerical  solution  of  Eq.  (2),  and  the 
particle  motion,  Eq.  (1),  establishes  the  particle  core  (PC)  model,  see  e.g.  [3].  As  we  are 
mainly  interested  in  the  influence  of  IBS  on  the  small  fraction  of  beam  particles  which 
may  form  a halo,  we  separate  the  collective  selfconsistent  space  charge  field  from  the 
individual  particle-particle  collisions.  The  basic  idea  is  now  to  use  the  PC  model  for 
a KV-beam  of  density  n and  charge  density  qii  as  a starting  point  and  to  superimpose 
a MD  treatment  for  a relatively  small  number  (~  some  hundred)  N,  of  test-particles 
in  a small  slice  moving  with  the  beam.  Applying  periodic  boundary  conditions  along 
the  beam  extension,  i.e.  in  z-direction,  the  test-particles  then  constitute  an  additional 
beam  of  density  ii,  -C  n which  is  subject  to  the  external  forces  K,  the  space  charge 
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fields  generated  by  the  underlying  KV-beam  and  the  direct  Coulomb  collisions  between 
the  test-particles.  But  in  the  test-particle  beam  IBS  is  strongly  reduced  as  compared  to 
the  real  beam  with  the  much  higher  density  n.  Because  no  length  scale  exists  for  the 
Coulomb  interaction,  we  compensate  the  smaller  test-particle  density  by  increasing  the 
test-particle  charge  to  a renormalized  value  The  N,  test-particles  with  mass  m and 
charge  qi  are  now  propagated  according  to  Eq.  (1)  feeling  the  same  external  forces  as 
the  physical  beam  particles  and  the  potential 


<5(ra,0 


N, 

y . 

p%4K£olr„-rp| 


J 4jceo|r„-r'| 


yab 


, (3) 


with  n = I / q^cnaobo,  oq  = a(s  = 0),  bo  = b(s  = 0)  and  the  step  function  0(z).  Here 
the  test-particles  directly  interact  as  charges  q,  but  couple  with  the  original  charge  q 
to  the  space  charge  field  of  the  underlying  real  beam.  To  avoid  double  counting,  the 
space  charge  field  (°c  q^n)  of  the  core  beam  is  corrected  by  the  test-particle  contribution 
(°c  q}n,).  This  combination  of  MD  simulations  and  the  PC  model,  the  particle-core- 
molecular-dynamics  (PCMD)  method,  allows  to  study  IBS  also  for  intensive  beams  on 
a qualitative  and,  at  least,  semi-quantitative  level,  with  a moderate  computational  effort. 
The  crucial  point  is  the  appropriate  choice  of  the  test-particle  charge  qt.  Taking  qi  —0 
we  retrieve  a pure  PC  description  with  a completely  two  dimensional  dynamics  in  the 
transversal  plane.  It  turns  into  a three  dimensional,  pure  MD  simulation  of  the  test- 
particle  beam  if  qfrit  approaches  q^n  and  the  space  charge  term  completely  vanishes. 
But  then  we  deal  with  a dilute  beam  of  high  charges  with  a much  lower  current  7,  = 
I{q,ni/qn).  The  disappearance  of  the  space  charge  term  »c  [q^n  - qjn,)  defines,  for  a 
given  density  n,,  a critical,  upper  limit  for  the  test-particle  charge  by  qc  = q{n/nty^^. 
Meaningful  applications  of  the  PCMD  description  are  restricted  to  values  q,  -C  qc,  where 
the  model  is  expected  to  reproduce  the  basic  effects  of  interparticle  correlations  and  IBS. 

The  coupling  within  a real  beam  of  density  n and  temperature  T is  described  by 
the  plasma  parameter  T = (g^/47teo^BT’)(4Jt«/3)'^^.  which  can  be  interpreted  as  the 
ratio  between  the  potential  energy  of  a pair  of  particles  at  their  mean  distance  and  the 
mean  thermal  energy.  For  a large  plasma  parameter  intrabeam  scattering  (IBS)  leads  to 
a secular  growth  of  emittance  [6].  But  here,  we  consider  beams  of  protons,  q — e,  with 
a large  temperature,  so  that  the  plasma  parameter  is  small,  with  typical  values  around 
r ~ 10~^.  It  is,  however,  a key  quantity  to  characterize  the  beam  properties,  and  the  test- 
particle  beam  in  the  PCMD  description  should  have  a similar  plasma  parameter  T,  as  the 
real  beam,  where  r,/F  = {qtlqfint/nyi^  = {qt/qc)^{qclqy^^  with  qc  = q{n/n,yi'^. 
An  analogous  argument  concerns  the  scattering  rate  v = «(v)o((v))  for  beam  particles 
with  an  averaged  relative  velocity  (v)  and  a scattering  cross  section  ct((v)).  It  can 
be  estimated  by  taking  the  Coulomb  cross  section  a ~ (<7^/47t£o«i(v)^)^.  This  yields 
V oc  nq'^  and  a ratio  of  the  scattering  rates  in  the  test-particle  beam  and  the  real  beam 
v,/v  = {qtlq^iit/n  = {qt/qcYiqdqy-  For  the  envisaged  high  current  proton  beams, 
test-particle  numbers  77,  ~ 1000  and  typical  settings  in  the  PCMD  simulations,  we  have 
7!,  ~ 10”^7rand^c~  1 0^  Hence,  a test-particle  charge  9,  ~ 10^  ~0.1^c  can  be  chosen 
which  allows  a matching  of  the  physical  relevant  parameters,  i.e.  F,  w F and  v,  a;  v, 
while  being  fairly  below  the  critical  value  qc  where  the  model  breaks  down. 
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PCMD  SIMULATIONS  OF  A BEAM  IN  A FODO  CHANNEL 


With  the  outlined  PCMD  scheme  we  studied  a continuous  beam  of  100  MeV  protons 
in  a FODO  channel  of  500  periods  of  1 m and  with  a beam  current  l=0.8A,  which 
may  represent  the  peak  current  in  the  bunch  for  a linac  with  a beam  current  in  the 
mA  range.  The  geometry  of  the  nominal  FODO  cell  was  taken  from  Ref.  [3].  The 
external  focusing  is  Kx  — -Ky  = 12m~-,  the  emittances  are  = e,.  = 10"^  m and  the 
space  charge  parameter  is  ^ = (e/)/(iteomc^y'P^)  = 9.6  x 10“^.  With  these  settings 
a matched  solution  of  the  envelope  equations  Eq.  (2)  is  constructed  numerically  and 
yields  «o  = 1.43  x 10“^  m,  = 8.54  x 10“'*  m,  dci/ds{0)  = 0 and  dh/ds{0)  = 0 at 
5 = 0 (t  = 0).  With  the.se  initial  condition.s,  the  envelopes  a{s),b{s)  are  propagated 
according  to  Eq.  (2)  throughout  the  structure  and  thereby  define  the  time  dependent 
core  field  contribution  in  expression  (3).  For  the  KV-distribution  the  corresponding 
beam  density  is  n = 1.0  x lO'^'cm'"^  and  the  plasma  parameter  is  F w 1.4  x 10 
when  averaging  over  the  temperatures  r,(.v)  and  Ty{s)  which  are  in  the  range  2.9  x 
10^ . . .8.4  X 10^  K and  are  connected  to  the  envelopes  via  ksTx/m  — e^(7Pc)^/4rt^  and 
ksTy/m  — The  initial  conditions  {rc((0),  v„(0)}  of  Nt  = 800  test-particles 

are  generated  according  to  a KV-distribution  in  the  transversal  degrees  of  freedom 
/(jc,y, v.v, Vy)  « b{\  — x^/a~  — y~/b^  — im’l/AkitTx  — inVy/^^ksTy)  and  homogeneously 
in  z-direction.  The  density  of  the  test-particle  beam  is  n,  = 7.3  x 10'*cm“-^  and  r/,.  = 
3.7  X 10^^.  The  simultaneous  propagation  of  the  V,  test-particles,  given  by  Eqs.  (1) 
and  (3),  is  treated  within  a usual  MD  scheme  by  standard  numerical  algorithms  with 
a time  step  which  guarantees  a sufficient  temporal  resolution  as  required  for  both  the 
rapid  changes  of  the  external  focusing  and  clo.se  Coulomb  collisions  between  the  test- 
charges.  In  z-direction  periodically  continued  slices  of  the  beam  are  considered,  while 
no  boundaries  exist  in  the  radial  direction.  To  study  a mismatched  beam,  a 20%  larger  «o 
than  for  the  matched  solution  is  chosen  and  taken  as  input  for  the  envelope  propagation 
(2)  with  the  same  space  charge  parameter  and  emittances.  Since  the  current  is  kept 
at  / = 0.8  A,  the  beam  density  n and  the  te,st-particle  density  n,  are  lower  for  the 
mismatched  beam,  but  the  ratio  n,/n  and  thus  are  the  same  as  for  the  matched  beam. 

In  our  first  studies,  we  concentrate  on  the  qualitative  effects  of  IBS  on  the  halo 
formation  and  are  mainly  interested  in  density  profiles  p(K,/)  and  particle  numbers 
N{K,t)  which  are  sampled  by  counting  and  normalizing  the  number  of  events  in  a given 
interval  of  the  normalized  transversal  position  K = Usually 

also  some  time  average  over  a couple  of  time  steps  is  performed. 

Unfortunately,  there  is  no  clear  cut  definition  of  a halo.  The  sharp  edges  of  a KV- 
beam  provide,  however,  at  least  a explicit  distinction  between  particles  inside  and  outside 
the  core.  This  is  already  anticipated  in  the  normalized  distance  K,  where  values  K < 1 
indicate  particles  inside  the  core,  regardle.ss  of  its  actual  shape  as  given  by  a{t)  and  /;(/). 
Fig.  1 presents  PCMD-simulation  results  of  the  particle  density  p(K,t)  as  a function 
of  K for  various  q,  and  different  positions  x in  the  FODO  channel,  and  for  a matched 
and  a mismatched  beam.  Clearly  visible  are  the  following  general  trends  of  the  density 
at  K > 1:  it  increa.ses  with  q,,  that  is,  the  strength  of  the  Coulomb  collisions,  i.e.  IBS, 
and  the  number  of  passed  structures,  and,  compared  to  the  matched  beam,  substantial 
larger  p(K  > 1)  are  observed  for  the  mismatched  beam.  These  effects  are  specifically 
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pronounced  at  the  largest  shown  = 50,  where  IBS  is  presumably  unrealistic  high.  The 
more  realistic  scenario,  with  T,  fa  F and  V;  ?a  v is  close  to  qt  = 10.  Here,  only  moderate 
effects  show  up.  The  difference  between  matched  and  mismatched  beam  is,  however,  as 
well  highly  significant.  The  relative  number  of  (test-)  particles  within  a certain 

interval  of  K > 1 (not  plotted)  shows  the  same  trends  as  p(k,t)  . Here  for  realistic  IBS 
rates,  corresponding  to  = 10. . .20  a fraction  of  ?a  lO""^  of  the  beam  particles  are 


0 0.5  1 1.5  2 0 0.5  1 1.5  2 

K K 

FIGURE  1.  PCMD  results  for  the  test-particle  density  p as  a function  of  the  dimensionless  distance 
K = {xl,fa{t)^  + ya/b{t)^y^^  for  various  test-particle  charges  q,  and  different  positions  t in  the  FODO 
channel  in  units  of  the  length  x,  of  one  elementary  structure.  The  matched  KV-beams  are  shown  on  the 
left,  the  mismatched  one  on  the  right. 
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around  k = 1 .5  at  the  end  of  the  channel  (t  > 400  tj).  But  we  remark  the  problem  of 
statistics  for  these  PCMD-runs  with  N,  = 800  particles  and  the  related  strong  fluctuations 
of  the  extracted  quantities.  When  performing  more  independent  simulation  runs  and/or 
using  larger  particle  numbers  we  expect  a comparable  fraction  of  particles  at  even  larger 
K. 

In  summary,  first  results  obtained  with  the  PCMD-model  for  a KV-beam  in  a FODO 
channel  demonstrate  that  even  rather  small  IBS  rates  significantly  enhances  the  number 
of  particles  outside  the  beam  core  in  a mismatched  high  current  beam.  In  the  matched 
beam,  yet  much  higher  IBS,  produces  only  a small  number  of  particles  slightly  above  the 
surface  of  the  core.  This  strongly  supports  the  existence  of  a scenario  whereby  already  a 
small  amount  of  IBS  is  sufficient  to  transport  a non-negligible  fraction  of  particles  in  the 
region  of  nonlinear  space  charge  forces  outside  the  beam  core  within  the  typical  time 
scale  of  beam  acceleration  in  a linac.  In  the  ca.se  of  a mismatched  beam,  resonances 
between  these  nonlinear  space  charge  forces  and  the  betatron  oscillations  push  the 
particles  away  from  the  core  to  higher  amplitudes.  This  mechanism  has  been  already 
identified  in  a pure  PC-picture  [3]  when  studying  particles  which  have  been  explicitely 
put  in  the  nonlinear  region  outside  the  core.  Including  IBS  the  population  of  this  region  is 
provided  by  the  hard  collisions  which  thus  serve  as  a trigger  for  the  elevation  of  particles 
by  the  space  charge  force  resonances.  This  supplies  a strong  evidence  that  IBS  serves  as 
an  indirect  source  for  the  halo  formation  in  a mismatched  beam.  To  verify  if  this  has  real 
significance  for  the  beam  losses  in  a high  current  linac,  subsequent  and  supplementary 
more  quantitative  investigations  are  needed. 
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